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A refined cascade model for kinetic turbulence in weakly collisional astrophysical plasmas is presented that 
includes both the transition between weak and strong turbulence and the effect of nonlocal interactions on the 
nonlinear transfer of energy. The model describes the transition between weak and strong MHD turbulence 
and the complementary transition from strong kinetic Alfven wave (KAW) turbulence to weak dissipating 
KAW turbulence, a new regime of weak turbulence in which the effects of shearing by large scale motions and 
kinetic dissipation play an important role. The inclusion of the effect of nonlocal motions on the nonlinear 
energy cascade rate in the dissipation range, specifically the shearing by large-scale motions, is proposed 
to explain the nearly power-law energy spectra observed in the dissipation range of both kinetic numerical 
simulations and solar wind observations. 



I. INTRODUCTION 

Plasma turbulence plays an important role in the 
transfer of energy in a wide range of space and astrophys- 
ical environments, from clusters of galaxies, to accretion 
disks around black holes, to the magnetized corona of our 
own sun, to the solar wind that fills the heliosphere. In all 
of these diverse environments, turbulence is responsible 
for the transport of energy from the large scales at which 
the turbulence is driven down to the small scales at which 
dissipation mechanisms can effectively dissipate the tur- 
bulent fluctuations and lead to heating of the plasma 
species. The understanding of plasma turbulence, its dis- 
sipation, and the resulting plasma heating is therefore a 
key goal for the space physics and astrophysics commu- 
nities. 

Although the study of the inertial range of plasma 
turbulence — the range of scales over which the effects of 
driving and dissipation are negligible — dates back more 
than four decades^, study of the dissipation range has 
only recently become a focus of the astrophysics and he- 
liospheric physics communities. One of the major chal- 
lenges in the investigation of the dissipation mechanisms 
at work in space and astrophysical plasmas is the fact 
that, at the small scales on which these mechanisms op- 
erate, the dynamics is typically weakly collisional. Un- 
der these conditions, a kinetic description of the plasma 
dynamics is necessary 2 - - — , a much more complicated de- 
scription than the standard fluid descriptions, such as 
magnetohydrodynamics (MHD), commonly used in the 
study of the inertial range dynamics. 

The increasingly vigorous activity studying the dissipa- 
tion of plasma turbulence has been driven by several new 
developments: the development of new analytical models 
of kinetic turbulences^, advances in high-performance 
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computation that have made possible kinetic simulations 
of weakly collisional turbulence^, and, most recently, 
new observational studies of turbulence in the solar wind 
with sufficient time resolution to extend the satellite- 
frame frequency power spectra of magnetic and electric 
field fluctuations up to the Dopplcr-shifted frequencies 
associated with the electron Larmor radius^—. 

This paper presents the refinement of a turbulent cas- 
cade model for weakly collisional plasmas, originally de- 
vised by Howes et al. 2008 5 . At the small scales where 
the turbulent motions are dissipated, two new physical 
effects may play a significant role: the transition between 
weak and strong turbulence, and the nonlocal contribu- 
tion to the nonlinear energy transfer rate. The aim of 
this paper is to describe how these additional physical 
effects are incorporated into a refined model, denoted 
the weakened cascade model. Energy spectra predicted 
by the new model are compared in detail to the numeri- 
cal energy spectra from nonlinear gyrokinetic simulations 
of turbulence in the dissipation range, and the physical 
implications of the nonlocal contribution to the nonlinear 
energy transfer in the dissipation range are discussed in 
detail. 



II. MOTIVATION FOR A REFINED CASCADE MODEL 

A new and exciting challenge in the study of turbu- 
lence in the solar wind is to understand the recent obser- 
vations of nearly power-law magnetic energy spectra in 
the dissipation ranged up to scales corresponding to the 
electron Larmor radius^™. The physics describing the 
turbulent dynamics of the inertial range is expected to be 
self-similar, dominated by the conditions that are "local" 
in scale 1 ®. In the dissipation range, on the other hand, 
dissipative mechanisms necessarily play a non-negligible 
role, and so the usual assumption that the turbulent en- 
ergy transfer is dominated by local interactions must be 
called into question. 
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In anisotropic Alfvenic turbulence, nonlinear interac- 
tions typically strengthen as the cascade progresses to 
smaller scale until they reach and maintain a state of 
critically balanced, strong turbulence. But as turbulent 
fluctuation amplitudes are diminished by some physical 
dissipation mechanism, one must consider the possibility 
that the turbulence may eventually transition back to a 
state of weak turbulence. In addition to this transition 
from strong to weak turbulence, one must also account 
for the effect on the nonlinear energy transfer rate at 
small scales (where the local turbulent fluctuation am- 
plitudes are diminished by some dissipation mechanism) 
by the undamped, nonlocal turbulent motions at larger 
scales. In this section, we briefly review the relevant con- 
cepts from modern theories of anisotropic magnetized 
turbulence that provide the theoretical foundation for 
these arguments. Then we identify the shortcomings of 
the original cascade model 5 - and outline the necessary re- 
finements of the model to account for the transition to 
weak turbulence at dissipative scales and for the effect of 
nonlocal motions on the energy cascade rate. 



A. Review of Anisotropic Plasma Turbulence Theories 

Early theories of magnetized plasma turbulence pro- 
posed an isotropic cascade of turbulent energy from low 
wavenumbers (large scales) to large wavenumbers (small 
scales) 17 ' 18 . Measurements of magnetized plasma turbu- 
lence both in the laboratory-i 9 - - — and in the solar wind 22 , 
as well as early numerical simulations^ 3 -, demonstrated 
the existence of significant anisotropy with respect to the 
direction of the local mean magnetic field. In the con- 
text of incompressible MHD, the inclusion of anisotropy 
in the direction of the nonlinear turbulent energy transfer 
through wavevector space lead to a fundamental distinc- 
tion between the properties of weak MHD turbulence and 
strong MHD turbulenc o 24 ' 25 . As the shift in behavior be- 
tween strong and weak turbulence plays a role when the 
turbulent cascade suffers dissipation, a brief review of the 
relevant properties is provided here. 

Although the theory for weak turbulence in incom- 
pressible MHD plasmas was a focal point of controversy 
in the late 1990fM^r^-, a refined theory of weak MHD 
turbulence has emerged. This theory is based on a per- 
turbative treatment of weak, resonant three-wave inter- 
actions between Alfven waves propagating in opposite 
directions along the local mean magnetic field. The need 
to satisfy both the resonance conditions for frequency 
and wavevector— and the requirement that only counter- 
propagating Alfven waves interact nonlinearl y 17 i 18 leads 
to the prediction that there is no parallel cascade of en- 
ergy in wavenumber space (where the directions parallel 
and perpendicular are defined with respect to the local 
mean magnetic field). Therefore, energy is transferred 
via the turbulent cascade anisotropically in wavevec- 
tor space, only to higher perpendicular wavenumbers 
k±, while the characteristic parallel wavenumber ku re- 



mains constant. The one-dimensional magnetic energy 
spectrum Es{kj_), defined such that the total magnetic 
energy Eb = J dk±EB(k±), is predicted to scale as 
Es(k±) oc kj 2 . Recent reduced MHD numerical sim- 
ulations have confirmed this predicted scaling in the 
regime of weak turbulence 3 ^.. The nonlinear interactions 
strengthen as the perpendicular wavenumber increases, 
so as the turbulent cascade progresses to higher fc^, the 
dynamics eventually violate the assumption of weak non- 
linear interactions required for the application of pertur- 
bation theory. Therefore, the theory of weak MHD tur- 
bulence predicts that, if a sufficiently large inertial range 
exists, weak MHD turbulence will eventually transition 
to a state of strong MHD turbulenc o 24 ' 25 . 

When the turbulence becomes sufficiently strong, the 
perturbative expansion of the nonlinear term is no longer 
dominated by the three- wave interactions, and terms of 
all orders contribute. This leads to a broadening of 
the resonance, relaxing the strict constraints on the fre- 
quency and wavevector, making possible the cascade of 
energy in the parallel direction of wavevector spaced 5 .. 
The theory of strong turbulence in incompressible MHD 
plasmas is based on the conjecture of a critical bal- 
anc o 25 ' 33 between the linear timescale for Alfven waves 
and the nonlinear timescale of turbulent energy trans- 
fer. The one-dimensional magnetic energy spectrum is 
predicted to scale as Es(k±) oc k^ ^ 3 . The conjec- 
ture of critical balance predicts the development of a 
scale-dependent anisotropy given by fc|| oc k 2 / 3 ; there- 
fore, the one-dimensional parallel magnetic energy spec- 
trum, defined by Eb — J dk^Esik^), is predicted to 
scale as Es(k\\) oc fc^ 2 . The predictions of the scale- 
dependent anisotropy are supported by early numerical 
simulation :) 34 ' 35 and by more recent studies of the scaling 
of the the parallel spectrum in the solar wind 36 ' 37 . 

The lack of a parallel cascade in weak MHD turbulence 
and the scale dependent anisotropy in strong MHD tur- 
bulence have an important consequence: even when tur- 
bulence is driven isotropically at low wavenumber (large 
scale) with fco = fc^o = fc||oi f° r a sufficiently large inertial 
range (typical of most space and astrophysical plasmas 
of interest), the turbulent fluctuations at high wavenum- 
ber (small scale) become significantly anisotropic with 
fen <C k±. When the turbulent cascade reaches the per- 
pendicular scale of the ion Larmor radius k±pi ~ I, such 
turbulent fluctuations transition to a cascade of strong 
kinetic Alfven wave (KAW) turbulence, as has been pre- 
dicted theoreticall y 5 ' 7 ' 38 ' 39 and verified with nonlinear 
kinetic simulations 8 . Assuming again a critical balance 
between the linear and nonlinear timescales, in the ab- 
sence of dissipation, the kinetic Alfven wave cascade is 
predicted to have a one-dimensional magnetic energy 
spectrum that scales as Es(k±) oc kj_ and a scale- 
dependent anisotropy k\\ oc k 1 / 3 for the nonlinear trans- 
fer of energy in wavevector space 5 -' 7 -' 4 ^— . When colli- 
sionless dissipation via the Landau resonance is included, 
the original cascade model 5 - predicts that, if magnetome- 
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FIG. 1. Schematic diagram of the anisotropic transfer of 
energy through wavevector space (kx,ku). For turbulence 
driven isotropically at low wavenumber k±op, ~ fenoPi ~ 10~ 6 
(quarter circle in the lower left of the figure), nonlinear inter- 
actions lead to transfer of energy to higher wavenumbers. For 
weak MHD turbulence, the transfer is restricted to higher k±, 
with no cascade of energy to higher parallel wavenumbers ku . 
As the MHD turbulence becomes strong (at k±pt ~ 10 -3 ), 
resonance broadening allows a parallel cascade of energy that 
achieves a state of critical balance at all scales, with k\\ oc k 2 / 3 . 
Finally, as the strong MHD turbulence transitions to strong 
KAW turbulence (at k±pi ~ 1), the turbulence remains in 
critical balance but the scaling in this regime changes to 
fey a k 1 / 3 ■ 



ter noise floor is taken into account, the spectral index 
of the measured one-dimensional magnetic energy spec- 
trum could vary from a value of —7/3 for weak dissipa- 
tion up to approximately —4 for strong dissipation, with 
the strength of the dissipation depending on the plasma 
parameters /3, and Tj/T e . 

To highlight the different behaviors predicted for weak 
MHD, strong MHD, and strong KAW turbulence, a 
schematic of the energy transfer through wavevector 
space is presented in Figure [TJ Since the dynamics of 
an MHD plasma is axisymmetric about the direction 
of the local mean magnetic field, the description of en- 
ergy transfer through the three-dimensional wavevector 
space can be reduced to the two-dimensional plane of 
the parallel component of the wavevector fcii (the ax- 
ial coordinate parallel to the local mean magnetic field) 
and the perpendicular component k± (the radial coor- 
dinate perpendicular to the local mean magnetic field). 
The anisotropic MHD turbulence theories above sug- 
gest that the energy flows through wavenumber space 
along the path given by the solid line in Figure [TJ — 
along constant hi for weak MHD turbulence and along 
a path given by critical balance for strong MHD and 
strong KAW turbulence. Numerical simulations of MHD 
turbulenc o 34 ' 35 ' 44 " — , however, suggest that the turbulent 
energy does not flow strictly along this one-dimensional 
path through wavevector space but approximately fills 
the shaded region in Figure [TJ Similar to the original 



cascade model, the one-dimensional weakened cascade 
model presented in this paper determines the turbulent 
energy integrated vertically in this figure over all possible 
values of kn at each k±; the effective value of kn for the 
integrated turbulent energy at each k± is assumed to be 
given by the solid line in Figure [TJ 

Several final comments regarding plasma turbulence 
theories are in order. First, note that the plot of the 
(k±,kn) plane in Figure [JJ is logarithmic on both axes. 
The inherently anisotropic transfer of energy in mag- 
netized plasma turbulence naturally leads to a condi- 
tion in which k\\ <C k± at high wavenumbers. A num- 
ber of studies of solar wind turbulence in the past have 
assumed a constant angle between the direction of the 
wavevector and the local mean magnetic field, which 
would correspond to a line with k\\ oc k± on the log- 
arithmic plot of wavevector space; this leads to a sig- 
nificant underestimate of the anisotropy at small scales 
(multi-spacecraft measurements in the solar wind demon- 
strate significant anisotropy at the scale of the ion Lar- 
mor radius^ 4 -). Second, although much of the weak and 
strong MHD turbulence theory applies rigorously only 
to incompressible MHD plasmas, it is our view that the 
fundamental concepts derived for incompressible MHD 
plasmas offer useful guidance for the understanding of 
turbulence in plasmas allowing a wider range of physi- 
cal effects, e.g., compressibility, finite Larmor radius ef- 
fects, linear kinetic damping. Finally, although we choose 
the particular strong turbulence scaling given by Gol- 
dreich and Sridhar— as the basis for the original and 
weakened cascade models, an alternative theory for the 
scaling of strong MHD turbulence has been proposed by 
Boldyre\« 4 i; application of the weakened cascade model 
using the Boldyrev scaling will be discussed in a subse- 
quent paper. 



B. The Original Cascade Model 

The aim of the original cascade model 5 is to explain 
the observed magnetic energy spectrum in the solar wind 
using a minimal number of ingredients, namely, finite 
Larmor radius effects and kinetic damping via the Lan- 
dau resonance. This model employs a one-dimensional 
continuity equation for the magnetic energy spectrum 
in perpendicular wavenumber space, and is based on 
three assumptions: (a) the Kolmogorov hypothesis that 
the energy transfer is determined locally in wavenumber 
spaced; (b) that a state of critical balance exists between 
the linear and nonlinear timescales at all wavenumbers in 
the spectru m 25 ' 33 ; and (c) that the linear kinetic damp- 
ing rates determine the dissipation of the turbulent fluc- 
tuations even in the presence of the nonlinear cascade. 
Results of the model generally demonstrate the qualita- 
tive feature that, as the linear kinetic damping becomes 
strong, the spectrum begins an exponential fall off. This 
qualitative feature is observed neither in recent nonlin- 
ear kinetic simulations^ nor in recent observations of the 
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dissipation range of the solar wind 1 ^— . 

The failure of the original cascade model to reproduce 
the correct qualitative behavior when collisionless damp- 
ing becomes strong is due to the assumption of strong 
turbulence satisfying the critical balance condition at 
all scales. It seems relatively clear that, as collision- 
less damping reduces the turbulent amplitude at a given 
wavenumber to a level below that expected for a dissi- 
pationless cascade, the nonlinear turbulent interactions 
at that wavenumber must cease to be strong. One then 
expects a transition to a cascade with weak nonlinear 
interactions (in the sense of weak turbulence), although 
the collisionless damping and shearing by larger scale mo- 
tions will also play an important role, so it won't be a 
standard weak turbulence picture. 

The failure of the strong turbulence assumption at per- 
pendicular wavenumbers where the dissipation becomes 
significant is made clear by a simple physical argument 
applied to the original cascade model results. In Fig- 
ure 7 of the original cascade model paper—, the parallel 
wavenumber as a function of the perpendicular wavenum- 
ber, k\\(k±), is plotted for plasma parameters ft = 1 and 
Ti/T e = 1. At a value of k±pi ~ 8, the value of ku peaks, 
and then it begins to drop. Physically, this would mean 
that, as the turbulence continues to cascade to smaller 
perpendicular scales, its parallel scale length actually in- 
creases. This predicted behavior does not make physical 
sense. This point is discussed at length in paragraph [46] 
of that paper, with the forward looking conclusion that 
"nonlinear simulations are necessary to determine accu- 
rately the behavior of the turbulent cascade as the kinetic 
damping becomes significant." 



C. Weak Turbulence and Nonlocal Interactions 

The more likely physical result, when the turbulence is 
dissipatively weakened, is that the cascade to higher par- 
allel wavenumber is suppressed, so the parallel wavenum- 
ber remains constant. Weak turbulence theory predicts 
that the parallel cascade is suppressed in incompressible 
MHD plasmas^, but whether this also holds in the ki- 
netic Alfven wave regime is not known. We conjecture 
here that the parallel cascade is suppressed for weak ki- 
netic Alfven wave turbulence, as discussed in more detail 
in ^III Al In addition to this effect, the nonlinear transfer 
of energy in weak turbulence also requires many uncor- 
related "collisions" between counter-propagating Alfven 
wave packets, in contrast to the single collision required 
in strong turbulenc o 24 ' 25 . 

Therefore, the weakened cascade model proposed here 
describes the transition between strong and weak turbu- 
lence by incorporating two changes: (1) altering the cas- 
cade in parallel wavenumber; and (2) increasing the num- 
ber of Alfven wave collisions required to accomplish non- 
linear energy transfer. The consequence of these changes, 
however, is that the collisionless damping at a given scale 
becomes relatively stronger than the nonlinear transfer, 



causing the spectra to cut-off even more abruptly, re- 
sulting in even greater discrepancy between the cascade 
model predictions and both the numerical simulations 
and the solar wind observations. 

To resolve this discrepancy requires the incorporation 
of another physical effect, in addition to the transition 
to weak turbulence. It requires accounting for the effect 
on the nonlinear energy transfer rate at a given scale by 
turbulent motions at other scales. Consider specifically 
the nonlinear turbulent energy transfer at a particular 
perpendicular wavenumber fc_L*. Although the turbu- 
lent energy is still transferred locally in scale space — 
for example, energy transfer from wavenumber kj_* to 
2k±* — this nonlinear energy transfer is not only due to 
motions at the local wavenumber fcj_», but also due to 
nonlocal (in scale space) motions at smaller and larger 
wavenumbers. Both coherent shearing by motions at low 
wavenumbers k± < fc^* (large scales) and incoherent dif- 
fusion by motions at high wavenumbers fcx > fc^* (small 
scales) may contribute substantially to the total nonlin- 
ear energy transfer rate from k±* to 2fcj_*. Therefore, 
we must abandon the Kolmogorov hypothesis of locality 
and account for the effect of nonlocal interactions on the 
turbulent energy transfer. 

The physical effect due to shearing by large scale mo- 
tions is essentially the same as that seen in high magnetic 
Prandtl number dynamo simulations using MHD 4 —. The 
magnetic Prandtl number is defined by the ratio of vis- 
cosity over the magnetic diffusivity, Pr m — v/rj, and typi- 
cally Pr m 3> 1 for many astrophysical plasmas of interest. 
For a high magnetic Prandtl number plasma, magnetic 
energy can be supported on subviscous scales, and this 
magnetic energy is observed in simulations to cascade to 
ever smaller scales, until the resistive scale is reached, at 
which point the magnetic energy can be dissipated. This 
cascade occurs although the plasma cannot support any 
fluid motions at the subviscous scales of interest. The 
transfer of magnetic energy at these small, subviscous 
scales is accomplished by the shearing due to larger scale 
motions (at scales larger than the viscous cut-off). This 
energy transfer is therefore nonlocal in nature. 

Such an effect must come into play in a weakly col- 
lisional plasma when the turbulent cascade reaches a 
small enough perpendicular scale that collisionless damp- 
ing can diminish the amplitude of the local (in scale) 
turbulent fluctuations. Kolmogorov's locality assump- 
tion is well supported in the inertial range, because the 
energy cascade rate due to shearing by local fluctuations 
always dominates over the energy cascade rate due to 
shearing by larger-scale motions. But, in the dissipation 
range, where the local fluctuations begin to diminish in 
amplitude due to some dissipative mechanism, their dom- 
inance of the local energy cascade rate breaks down, and 
the effects of shearing by larger-scale, undamped motions 
must be taken into account. In this way, the energy cas- 
cade rate at the small, rather strongly damped scales may 
be dominated by this nonlocal shearing, with the weak 
turbulent interactions due to local fluctuations playing a 
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subdominant role. Thus, the energy can be cascaded to 
ever smaller scale at a reasonably large rate even though 
the nonlinear energy transfer rate due to the local fluctu- 
ations becomes negligible. This nonlocal shearing effect 
is proposed to explain the nearly power-law appearance 
of the numerical and observational spectra, and is critical 
for understanding the turbulence in the dissipation range. 
The nonlocal contribution to the energy cascade rate by 
diffusive motions at small scales, although negligible at 
the far end of the dissipation range, is also included in 
the refined model for consistency. 



III. THE WEAKENED CASCADE MODEL 

In order to refine the original cascade model*, we must 
incorporate two physical effects not included in the orig- 
inal model: 

1. Weak Turbulence: The model must be able to 
handle the quantitative changes in the energy cas- 
cade rate in both the weak and strong turbulence 
regimes. Of particular importance is the qualita- 
tive difference in the parallel cascade of energy: in 
weak turbulence, there is no cascade of energy to 
smaller parallel scales; and in strong turbulence, 
the parallel cascade is governed by critical balance. 

2. Nonlocal Interactions: The net energy cascade rate 
at a given wavenumber e(k±*) must account for the 
nonlinear transfer due to both the local fluctuations 
at kj_* and nonlocal fluctuations at other wavenum- 
bers. The nonlocal fluctuations contribute to the 
nonlinear energy cascade rate due to shearing by 
fluctuations at smaller wavenumbers fcj_ < 

and diffusion by fluctuations at larger wavenum- 
bers k± > 

In this section, we describe in detail the quantitative 
modifications of the original cascade model required to 
incorporate these two physical effects, resulting in the 
new weakened cascade model. 



A. Weak Turbulence 

The key parameter in distinguishing weak from strong 
turbulence is the nonlinearity parameter 



k±vj 



X 



(1) 



This dimensionless parameter x measures the ratio of the 
nonlinear frequency ui n i ~ kj_v± to the linear wave fre- 
quency ui. Strong turbulence corresponds to \ ~ 1, sat- 
isfying the condition of critical balance ui ~ ui n i , whereas 
weak turbulence corresponds to x C 1. Note that the 
case of overdriven turbulence \ > 1 has not been thor- 
oughly explored or discussed in the literature; it probably 
deserves some attention, but henceforth we will consider 



only the cases of weak or critically balanced turbulence 

To handle the transition to kinetic Alfven wave tur- 
bulence, we focus on the magnetic energy rather than 
the kinetic energy, so we want to write the nonlinear en- 
ergy cascade rate in terms of the magnetic fluctuation 
energy, b\ = 6B 2 L (k±)/4irriimi J where we have written 
the magnetic field fluctuation in velocity units. We also 
adopt the shorthand v k = v±(k±). Following the original 
cascade model, we relate the velocity and magnetic field 
fluctuations to each other using the linear theory, 



v k = ±a{k 1 _)b k 



(2) 



where the coefficient a smoothly transitions from the 
MHD to the KAW limit, 



a(kj_) 



I . k±pi < f 

k_LPi/y/Pi + 2/(l+T e /Ti), fc ±Pl »I 



(3) 



Note that the linear frequency in the gyrokinetic limit 49 
is given by 



(4) 



where uj(k±) — a(k±) in both asymptotic ranges k±pi <C 
I and k±pi 3> I but not in the transition region k±pi ~ 1. 
For this simple model, we take the approximation that 
a(k±) =uJ(k±) over all scales. 

Thus, an appropriate definition of the nonlinearity pa- 
rameter valid for all scales is 



C 2 k±v k C 2 k±b k a C 2 k±b k 



k\\VACt 



k\\VA 



(5) 



where C 2 is an order-unity, dimensionless "Kolmogorov" 
constant. 

There are two primary effects that must be captured 
in order to incorporate the weak turbulence limit into the 
cascade model: 

1. Slower nonlinear energy cascade rate due to the 
effect of many uncorrelated weak interactions be- 
tween oppositely directed Alfven waves^. 

2. The suppression of the energy cascade to small par- 
allel scales. 

These issues are addressed separately below. 
We define the energy cascade rate by 



e(fcjj = C\ 3/2 u) nl b 2 k . 



(G) 



where C\ is an order-unity, dimensionless Kolmogorov 
constant and the nonlinear frequency, bridging weak and 
strong turbulence, is given by 

u n i = xk±v k = ( -J—^h. ] k^b k a. (7) 
V k \\ v A J 
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Note that the cascade time is given by 

knV A 



1~nl 



C 2 k±b k 



(k±b k a) 1 



(8) 



k\\VA 



C 2 {k\\v A a)- 1 = N 2 C 2 (k llVA a)-\9) 

C 2 k\b k J 

This shows explicitly that, in the MHD limit a = 1, 
it takes N 2 = 1/x 2 Alfven wave packet collisions, each 
lasting an Alfven wave crossing time ta = (k^VA) -1 , for 
the energy at a given scale to be transferred to the next 
scale. 

In summary, the energy cascade rate may be written 

' C 2 kxb k 



_ n -3/2 



k\\v A 



k±bia. 



(10) 



Let us investigate the predicted steady-state magnetic 
energy spectrum in different limits for a constant energy 
cascade rate e. 

In the MHD limit, a — 1. The critically balanced, 
strong turbulence limit gives \ = C 2 k±bk/k\\VA = 1, 

so we have e — C\ 3 ^ 2 k±bf,, yielding a solution for 

the magnetic field b k = C\^ 2 'e 1 / 3 k A 1 ^ 3 and a corre- 
sponding 1-D magnetic energy spectrum £^B(fcj_) = 

bf./k± = Cie 2 / 3 k ± 5 ^ 3 . In the weak turbulence limit 

X < 1, we have e = C 1 3 ^ 2 C 2 k 2 1 b\jknVA and ku re- 
mains constant, so we obtain a solution for the mag- 
netic field b k = Cl /8 C 2 ~ 1/4! e 1 / i (k l \VA) 1/i kl' L/2 and a cor- 
responding 1-D magnetic energy spectrum Es(k±) = 

C^C-^e^mvA) 1 / 2 ^ 2 . 

In the KAW limit, we approximate a ~ k±pi. 
Strong, critically balanced KAW turbulence also has x = 
C 2 k±bk/k\\ va = 1, so we have e = C 1 3 ^ 2 'k\pib\, yielding 

a solution for the magnetic field bk = C\^ 2 t^l^k 2 p i 1 ^ 3 
and a corresponding 1-D magnetic energy spectrum 
E B (k ± ) = C 1 e 2 / 3 k~ 7/3 p~ 2/3 . In the weak turbulence 

limit x < 1) wc have e = C\ 3 ^ 2 C 2 k\b\pijk\\y a- In 
this case, we assume that fen remains constant in the 
weak turbulence limit; this assumption is discussed later 
in this subsection. The solution for the magnetic field 
is then b k = C 3/8 C 2 - 1/4 e 1 /4( fc||l;A )i/4 fc -3/4 p -i/4 and & 

corresponding 1-D magnetic energy spectrum EB{kj_) — 
C 3 C 2 1 ' 2 e 1 / 2 (k\\VA) 1 / 2 kJ 5 p^ 1 ^ 2 agreeing with theo- 
retical predictions for weak KAW turbulence^. 

The other important effect when synthesizing a theory 
that combines the weak and strong turbulence limits and 
their influence on the turbulent cascade of energy is to 
model the parallel cascade of energy. In the strong limit, 
the condition of critical balance allows the determination 
of k\\(k±), strictly a function of k±. The assumption of 
critical balance at all scales in the original cascade model 
allowed the damping term to be simplified. The original 
cascade model equation for the evolution of the magnetic 
energy was written as 

' ; ' ie(fc -il + 5(fc ± )-2 7 6 2 , (11) 



dt 



where the last term could be simplified, under the as- 
sumption of critical balance and for the linear kinetic 
damping rate in the gyrokinetic limit 7 = i^(kj_)k\\VA, 
to take the form 



w(/c_l) 



(12) 



strictly a function of k±. When this assumption is re- 
laxed, however, to determine the damping rate one needs 
to know the value of fc|| (fc_i_) to find 



7(fe_L,fc||) = j(k±)k\\v A . 



(13) 



Without critical balance, we must devise some other 
means of determining an appropriate value of fc|| . 

It must be noted that the distribution of turbu- 
lent power in the plasma can fill a region of the two- 
dimensional wave vector space (fcj_,A;||), as denoted by 
the shaded region in Figure [TJ Here we follow the orig- 
inal cascade model in treating the magnetic fluctuation 
energy integrated over all possible values of fc|| , so that 



b'i(k±)= I dk\\SB^(k±,k\\)/4imimi. (14) 



We will treat the total magnetic energy as if it resides at 
a single value of k\\ . The particular choice we make is the 
maximum value of fcii that contains significant fluctua- 
tion energy. Therefore, in the strong turbulence limit, k\\ 
is chosen so that the parallel cascade satisfies the critical 
balance criterion. In the weak turbulence limit, how- 
ever, we assume that the parallel cascade of energy is 
inhibited, so the value of ku remains constant as energy 
cascades to larger k±. By making this choice, we again 
can solve for the turbulent cascade of magnetic energy as 
a one-dimensional problem in k± , but as we shall see this 
necessitates solving for fcii from the driving wavenumber 
on up. 

It has been proven rigorously that there is no paral- 
lel cascade of energy for weak turbulence in the limit 
of incompressible MHD^. Heuristic arguments for weak 
turbulence in incompressible Hall MHD plasmas 50 and 
numerical evidence demonstrating an anisotropic cascade 
of energy in weak whistler wave turbulence^— suggest 
that the parallel cascade is also suppressed for the dis- 
persive wave modes at perpendicular scales smaller than 
the ion Larmor radius. Therefore, we conjecture that 
the parallel cascade is suppressed in weak kinetic Alfven 
wave turbulence as well^. Nonlinear kinetic simulations 
of turbulence in the kinetic Alfven wave regime will play 
a key role in testing this hypothesis. 

To model these effects, we take the equation for the 
evolution of the parallel wavenumber to be 



rfln/c 



d In fcj 



2/3 + (l/3)(fc^) 2 
1 + (k± Pl ) 2 



(15) 



dk x 



Let us now consider the limits of this equation. In the 
limit of critically balanced, strong turbulence, x = 1- In 
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the MHD limit kxPi <C 1, we obtain dhiku/dhikx = 
2/3, or fcii oc k 2 ^ 3 ; in the KAW limit, kxPi 1, we 
find dbxkn/dhikx = 1/3, or ku oc k 1 / 3 . In the weak 
turbulence limit in both regimes, \ 2 ~^ 0, so we find the 
result that k\\ remains constant^. 

To apply this equation to the cascade model, we begin 
at the driving scale (fcj_o, fc||o)) an d integrate forward over 
the logarithmically spaced grid points in k± 



k h+l = fc| U |l + 



2/3+(l/3)(fc ±jft ) 2 
1 + (k±jpi) 2 



(16) 

where the nonlinearity parameter must be calculated at 
each gridpoint j, 



X.i 



hjVA 



(17) 



Note that the value of the nonlinearity parameter is con- 
strained to have values Xj < 1 at all points. 

Note also that ku is never allowed to decrease as k± 
increases. This behavior appears to make good sense 
physically. Consider turbulent fluctuations at a given 
scale characterized by a perpendicular scale lx* (or kj_*) 
and a parallel scale hu (or fciu). For a fluctuation at 
a smaller perpendicular scale l± < lx* (or kx > kx*), 
it seems to be unphysical that this smaller perpendicular 
scale fluctuation could have a larger parallel scale In > in* 
(or fc|| < 



B. Nonlocal Interactions 

The Kolmogorov hypothesis states that the rate of non- 
linear energy transfer at a given perpendicular wavenum- 
ber kx* depends only on the conditions at that wavenum- 
ber. This assumption of locality in scale space leads to 
the familiar self-similar scaling of the turbulent cascade 
within the inertial range. When the effects of dissipation 
are taken into account, however, it becomes necessary 
to abandon this limiting hypothesis and to account for 
the effect of nonlocal motions, at both lower and higher 
wavenumbers, on the nonlinear energy cascade rate at 
the local wavenumber kx*- 

We first note that the energy transfer still occurs lo- 
cally in scale space, with energy being transferred from 
wavenumber kx* to 2k x*- This energy transfer, how- 
ever, is not mediated solely by the local turbulent mo- 
tions at wavenumber kx*- Rather, motions from both 
smaller wavenumbers k± < kx* (large scales) and larger 
wavenumbers kx > kx* (small scales) may contribute 
substantially to the total nonlinear energy transfer rate 
at wavenumber kx*. Therefore, the cascade described 
by the weakened cascade model is essentially nonlocal 
in character, formally precluding the possibility of self- 
similar solutions. Below we discuss the effect of both 
larger and smaller scale motions on the local turbulent 
energy transfer. 



A generic property of the turbulent cascade — whether 
weak or strong or in the MHD or KAW regimes — is that 
the nonlinear timescale t„j ~ w"; 1 always decreases as the 
wavenumber increases. Therefore, the longer characteris- 
tic timescale of large scale turbulent motions (relative to 
the local scale) means that their effect on the local tur- 
bulent fluctuations will be coherent over the lifetime of 
the local fluctuations. The large scale turbulent fluctua- 
tions can effectively be considered as a shearing motion 
applied to the local fluctuations. The contribution to 
the nonlinear frequency due to these large scale shearing 
motions can be accounted for by summing over all larger 
scale motions in the cascade, 



» 



(fcj 



d]nk' ± J? C) (k' ± )e(kx-k' ± ), (18) 



where kxo is the driving (outer) scale of the turbulence, 
k±max is the smallest (inner) scale, and is the piecewise 
constant Heaviside step function. Here, the nonlinear 
frequency due to shearing motions at a scale k'j_ is given 
by 

^° C) (fcl) = X(k' ± )k' ± b k (k' ± )a(k' ± ). (19) 

The effect of smaller scale motions on the local scale, on 
the other hand, can be treated as a diffusive process be- 
cause the characteristic lifetime of the small scale motions 
is shorter than the timescale of the local fluctuations. 
The diffusion coefficient due to motions at a scale I with 
a timescale t\ is given by T> ~ Z 2 /r;. When translated 
into nonlinear frequencies and wavenumbers, the diffu- 
sion coefficient due to motions at a large wavenumber 



k' ± (small scale) takes the form T>(k' ± ) 



(loc) 

J nl 



(k' ± )/k 
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Summing over all smaller scales leads to the contribu- 
tion to the nonlinear frequency due to these small scale 
diffusive motions, given by 

o#(fc ± ) = f k ' maX dlnk^i^ik'^Qik^ kx). 

( 20 ) 

Therefore, the total nonlinear frequency may be writ- 
ten as the sum of terms due to large-scale shearing, local- 
scale fluctuations, and small-scale diffusion, 



u> n i{kx) = ^(kx) + + (21) 

where the local term may be written in an analogous 
manner, 



,(') 



d\nk' ± J^ C) (k' ± )S(kx - k'x). (22) 



It is important to note that the inclusion of these non- 
local interaction terms does not alter the scaling of the 
turbulent cascade in the MHD inertial range (or in the 
KAW inertial range in the absence of dissipation), as 
shown in Appendix [Al 
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C. Summary of the Weakened Cascade Model 

In this section, we summarize the weakened cascade 
model for clarity and ease of reference. The key depen- 
dent variables of this one-dimensional model are the per- 
pendicular magnetic field energy b\ (integrated over all 
feii) and the parallel wavenumber fen, both functions of 
the independent variable, the perpendicular wavenum- 
ber The continuity equation for magnetic energy in 
perpendicular wavenumber space is given by 



dt 



(23) 



where the steady state magnetic energy spectrum is de- 
termined by iterating numerically until the right-hand 
side equals zero for all wavenumbers. Here the source 
term S(k±) determines the energy input at the driving 
scale, characterized by driving wavenumber components 
fc^o and fen q . The linear kinetic damping rate is a func- 
tion of both fc_L and ku and may be written as 



7(fej_,fe||) = j(kx)k\\v A - 



(24) 



Note that, in the gyrokinetic limit fe|| <C fcj_, the nor- 
malized damping rate 7 = j/(k\\VA) is only a function of 
fej_ and may be determined by the gyrokinetic or Vlasov- 
Maxwell dispersion relation 5 -^. The energy cascade rate 
e(fe±) is a subsidiary function defined by 



e(fcj_) = C 1 3/2 0J n ib 



(25) 



and the total nonlinear frequency uj n \ at wavenumber 
k± — including terms due to large-scale shearing, local- 
scale fluctuations, and small-scale diffusion — is given by 

f k _|_ m a x 



w n /(fcj_) 



d In k'j_t 



floe) 



'nl 



e(fcj 



k' ± ) + ^Q(k' ± -k ± ) 



(26) 



_i_ —i 

where the piecewise constant Heaviside step function is 
defined by 



Q(x) 



1 x > 
1/2 x = 
x < 



(27) 



and where the contribution to the nonlinear frequency 
due to motions at each wavenumber k' ± is given by 



^ c \k' 1 _) = x (k' 1 _)k'Mk' 1 _p(k' ± ) 1 

and the nonlinearity parameter is defined by 

C 2 k^b k {k^ 



X(k'± 



min 1 



fc||(fe^_)i;^ 



(28) 



(29) 



In the gyrokinetic limit feii <C k±, the linear Alfven wave 
frequency is given by uJ(kj_) — u / (k\\VA) and may be de- 
termined by the gyrokinetic or Vlasov-Maxwell disper- 
sion relation 5 -! 5 - 6 -. Finally, the value of fci|(fc_i_) is found 



by integrating from the driving scale (fe^O; fe||o) using the 
equation 



din fen 
din k± 



2/3+(l/3)(fc_LPi 



1 + (k±p*) 2 



(30) 



The model has only two free parameters, the dimension- 
less order-unity Kolmogorov constants: C\ adjusts the 
relative weight of the nonlinear energy transfer to the 
linear kinetic damping, and C2 adjusts the condition of 
critical balance. 



IV. NUMERICAL SIMULATIONS 

To test the ability of the weakened cascade model to 
predict the properties of the turbulent steady-state in 
a weakly collisional plasma, we employ nonlinear gyroki- 
netic numerical simulations using the code AstroGK^. In 
particular, we are interested in the ability of the weak- 
ened cascade model to reproduce the turbulent energy 
spectra when the dissipation becomes strong, so we focus 
on two simulations that each cover the entire dissipation 
range of scales from the ion to the electron Larmor ra- 
dius. For one simulation we choose parameters /?; = 1 
and Ti/T e = 1, leading to moderate collisionless damp- 
ing of the turbulent cascade; for the second simulation we 
choose parameters /3; = 0.01 and Ti/T e — 1, leading to 
strong collisionless damping. All the details of the — 1 
simulation are outlined in Howes et al. 2011—, so here we 
focus on the description of the new = 0.01 simulation. 

AstroGK evolves the perturbed gyroaveraged distribu- 
tion function h s (x, y, z, A, e) for each species s, the scalar 
potential tp, parallel vector potential Au, and the par- 
allel magnetic field perturbation <5-B|| according to the 
gyrokinetic equation and the gyroaveraged Maxwell's 
equation s 56 ' 58 . The velocity space coordinates are A = 
v\/v 2 and e — v 2 /2. The domain is a periodic box of 
size L\ x L|| , elongated along the straight, uniform mean 
magnetic field Bq. Note that, in the gyrokinetic formal- 
ism, all quantities may be rescaled to any parallel dimen- 
sion satisfying Ly/L^ 3> 1. Uniform Maxwellian equilib- 
ria for ions (protons) and electrons are chosen, and the 
correct mass ratio mi/m e — 1836 is used. Spatial dimen- 
sions (x,y) perpendicular to the mean field are treated 
pseudospectrally; an upwinded finite-difference scheme 
is used in the parallel direction, z. Collisions are in- 
corporated using a fully conservative, linearized collision 
operator that includes energy diffusion and pitch-angle 
scatterin g 59 ' 60 . 

Both simulations employ a simulation domain size 
L±_ = 2npi with dimensions (n x ,n y ,n Zl n\,n £ ,n s ) — 
(128, 128, 128, 64, 16, 2). The fully dealiased range yields 
perpendicular wavenumbers 1 < k x pi < 42 and 1 < 
kyPi < 42. For a mass ratio of m,i/m e — 1836 and 
Ti/T e = 1, the scale of the electron Larmor radius fcj_p e = 
1 corresponds to a value of k±pi ~ 42.8, so the simulation 
covers the entire dissipation range of scales from the ion 
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to the electron Larmor radius. The importance of cov- 
ering this entire range is that wave-particle interactions 
via the Landau resonance are resolved and sufficient to 
damp the electromagnetic fluctuations within the simu- 
lated range of scales. Therefore, all dissipation in the 
simulation is accomplished by resolved Landau damping, 
eliminating the need for an ad hoc fluid model of the 
damping, such as viscosity or resistivity 6 . 

Both simulations are driven at only the largest per- 
pendicular scale in the domain, k±pi = 1, using six 
modes of a parallel "antenna" current j| k added via 
Ampere's Law££. For the Pi = 0.01 simulation, these 
driven modes have wavevectors (k x pi,k y pi,k z L\\/27r) = 
(1, 0, ±1), (0, 1, ±1), (-1, 0, ±1), frequencies uj a = 1.2w A0 
(where u>ao = k\\oVA is a characteristic Alfven frequency 
corresponding to the parallel size Lu of the domain), and 
amplitudes that evolve according to a Langevin equation. 
This produces Alfvenic wave modes with a frequency 
lu ~ ±fc|| ?M and a decorrelation rate comparable to to, 
as expected for critically balanced Alfvenic turbulence 25 . 

The coefficients for the collision operato r 59 i 60 in the 
Pi = 0.01 simulation are vi — O.Olu^o and v e = 
O.loMoj chosen to achieve sufficient damping of small- 
scale velocity-space structure yet to avoid altering the 
collisionless dynamics of each species over the range of 
scales at which the kinetic damping is non-negligible. 

Both simulations are brought to a statistically steady 
state with minimal computational expense by using a re- 
cursive expansion procedure^. At low spatial resolution, 
the simulation is run for more than an outer-scale eddy 
turnover time (this turnover time is To = 4.65oj^q for the 
Pi = 0.01 simulation) to reach a steady state. Resolu- 
tion in each spatial dimension is then doubled, and the 
simulation is run to a new steady state, which requires 
only a time of order the cascade time at the smallest re- 
solved mode before expansion. Three applications of this 
expansion procedure brings the perpendicular dynamic 
range from an initial value of 6 to a final value of 42. 
The Pi = 0.01 simulation has been evolved using this 
procedure to a time of t = 166. 230^. 

The resulting magnetic and electric energy spectra 
from both simulations are presented in the next sec- 
tion (Figures [3] and [5]) . The normalized one-dimensional 
magnetic-energy spectrum is defined by EB ± (k±) = 
(L ll /L 2 L )2n(k ±Pl ) 3 J dz (\A lk± (z)\ 2 ) /Snn.T,, where the 
angle brackets denote a sum of the energy of all perpen- 
dicular Fourier modes falling in a wavenumber shell cen- 
tered at |kj_| = k± with width 2t:/Lj_. The normalized 
electric-energy spectrum EE ± (kj_) is defined similarly in 
terms of ^k ± , with an extra factor of (c/va) 2 , where c 
is the speed of light. Note that since the Fourier modes 
are defined in Cartesian coordinates, but the value of 
kj_ is polar, there are modes in the "corner" of the sim- 
ulation (in Fourier space) where full information at all 
azimuthal angles is not available — therefore beyond the 
value of k±pi = 42, not all modes are represented, lead- 
ing to a drop in the energy in the 1-D energy spectrum. 
We choose to plot the spectra over this range to demon- 



strate that no bottleneck of electromagnetic fluctuation 
energy at small scales occurs in these simulations. 



V. RESULTS OF WEAKENED CASCADE MODEL 

In this section we present results of the weakened cas- 
cade model, as defined by equations (j2"3"]) - (j3"0)) . for cases 
of interest. All weakened cascade model calculations pre- 
sented in this section employ the linear collisionless gy- 
rokinetic dispersion relation^ to calculate the normalized 
linear kinetic frequency ZJ = uj /(k\\VA) and damping rate 
7 = /(k\\VA) as a function of three dimensionless plasma 
parameters: the normalized perpendicular wavenumber 
k±pi, the ion plasma beta Pi, and the ion to electron 
temperature ratio Ti/T e . A fully ionized plasma of pro- 
tons and electrons with isotropic Maxellian equilibrium 
velocity distributions is assumed, and a realistic mass ra- 
tio of rrii/m e = 1836 is used. 



A. Transition from Weak to Strong MHD Turbulence 

The first important test of the weakened cascade model 
is to verify that it reproduces the theoretically predicted 
characteristics of the transition from weak MHD turbu- 
lence to strong MHD turbulence. The plasma parameters 
chosen for this test are Pi = 1 and Ti/T e = 1, and the tur- 
bulence is driven isotropically at k±opi — k^ pi = 10~ 6 . 
The amplitude of the source term is chosen so that the 
nonlinearity parameter at the driving scale is x(^J-o) = 
0.1, giving rise to a weak MHD turbulent cascade. The 
Kolmogorov constants are taken to be C\ = 1.4 and 
C*2 = 1, and the spectrum is solved using 60 logarith- 
mically spaced gridpoints over k±pi G [10~ 6 , 1]. 

The results of this test, presented in Figure [H are con- 
sistent with theoretical predictions, showing a transition 
from the characteristics of weak MHD turbulence at low 
perpendicular wavenumber to those of strong MHD tur- 
bulence at high perpendicular wavenumber. Presented 
in the panels of Figure [5] are: (a) the energy spec- 
trum of the perpendicular magnetic field fluctuations, 
Eb ± = b 2 /k±, normalized by the value of the spec- 
trum at the driving scale, (b) the normalized parallel 
wavenumber fcyp^, and (c) the nonlinearity parameter \i 
all plotted vs. the normalized perpendicular wavenum- 
ber k±pi. In the weak MHD turbulence regime, over the 
range of scales 10~ 6 < k±pi < 10~ 4 , the magnetic en- 
ergy spectrum yields the theoretically predicted scaling 
Eb ± oc kj 2 , the parallel wavenumber remains essentially 
constant, and the value of the nonlinearity parameter 
is small with x < 0.5. In the strong MHD turbulence 
regime 10~ 3 < k±pi < 1, the magnetic energy spectrum 

— 5/3 

scales as Eb ± oc k ± , the parallel wavenumber scales 

as predicted by the critical balance hypothesis kit oc k 2 ^ 3 , 
and the nonlinearity parameter reaches and maintains a 
value x ~ 1- Note that the nonlinearity parameter is 
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FIG. 2. Results of the weakened cascade model for the tran- 
sition from weak MHD turbulence to strong MHD turbulence 
for a plasma with f3i = 1 and Ti/T e = 1. Plotted in the panels 
are: (a) the logarithm of the energy spectrum of the perpen- 
dicular magnetic field fluctuations, log(EB ± ) = log(&|/fcj_), 
normalized by the value of the spectrum at the driving 
scale, (b) the normalized parallel wavenumber k^pi, and (c) 
the nonlinearity parameter \ vs - the normalized perpendic- 
ular wavenumber k±pi. Results show a weak MHD turbu- 
lence regime for k±pi £ [10 -6 ,10~ 4 ], a transition regime for 
k±pi € [10 -4 , 10 -3 ], and a strong MHD turbulence regime for 
k±pi £ [10 -3 , 1]. 



numerically constrained to be x < 1 (dashed line); cal- 
culating x without this constraint yields the dotted line 
shown in panel (c), which still remains of order unity, 
with x 1-5- Each of these results for the weak and 
strong MHD turbulence regimes agree with theoretical 
predictions. The range of scales 10~ 4 < k±pi < 10~ 3 is 
a transition region that smoothly connects the weak and 
strong MHD turbulent cascades. 

Given that the state of the art in numerical simulations 
of turbulence can reach only a factor of approximately 
10 3 in each dimension, these results would suggest that 
it is not currently possible to capture the physics of the 
transition from the weak to the strong MHD turbulence 
regime in a single numerical simulation. Since the transi- 
tion regime spreads over a factor of 10 in dynamic range, 
and one would need a minimum dynamic range of 10 
in each of the weak and the strong turbulence regimes, 
these ranges alone would require all of the resolution cur- 
rently feasible, leaving no room for energy injection and 
dissipation. Recent work simulating sub-ranges of this 
plot, however, have indeed confirmed the transition from 
the weak to the strong scaling of the magnetic energy 
spectrum in a series of reduced MHD simulations 3 ^. 



B. Local vs. Nonlocal Models 

To test the importance of the nonlocal contribution to 
the turbulent energy cascade rate, we can define a com- 
parable local model by defining the nonlinear frequency 
at a given perpendicular wavenumber, 

Wnj(fcx) = C 3 J^ c) (k ± ). (31) 

Note that the constant C3 does not represent an addi- 
tional free parameter for this local model — this constant 
may be absorbed into the Kolmogorov constant C\ — but 
is introduced to enable easy comparison with the non- 
local model given by equation (126[) . The integration of 
equation (l2"fJl) in the strong MHD and strong KAW in- 
ertial ranges (when damping is negligible), as presented 
in Appendix [A] suggest that, to yield comparable local 
and nonlocal models, the value of this constant should 
be set to C3 = 2.25; in this case, the values of C\ and 
C2 in the local and nonlocal models should be directly 
comparable. We shall see below that, although the local 
model can fit nonlinear numerical simulation results in 
certain circumstances, it is missing the essential physics 
required to fit a wide range of cases. 

In evaluating the nonlocal model, it is desirable to de- 
termine the relative contributions of large-scale shear- 
ing, local scale motions, and small-scale diffusion to the 
nonlinear frequency at particular wavenumber, ui n i{k±). 
If the definition of the local contribution given in equa- 
tion (f2"2"j) is used, the delta function ensures that the local 
contribution is an infinitesimal slice of the entire integral, 
and the local contribution is not easily comparable to the 
large- or small-scale contributions. To yield more eas- 
ily interpretable results for the contributions to u> n i(k±), 
we split the integral in equation (|26l) into three ranges: 
the large-scale shearing contribution (s) over [fc_i_o, k±/2), 
the local-scale contribution (1) over [k±/2,2k±], and the 
small-scale diffusive contribution (d) over (2fcj_, k± rnax }. 

Although the weakened cascade model follows only 
the cascade of perpendicular magnetic energy, the en- 
ergy spectra of other fields can be constructed from the 
stready-state solution. Assuming the waves have the 
character of the linear Alfvenic eigenmodes, we use the 
solution of the linear kinetic eigenfunction as a function 
of k±_ to construct, for example, the amplitude of the per- 
pendicular electric field fluctuation from the amplitude of 
the perpendicular magnetic field fluctuation given by the 
cascade model solution. Because the phase and ampli- 
tude relations between the fields are fixed by the linear 
kinetic physics, no additional free parameters are intro- 
duced: if the linear character of the fluctuations applies, 
the solution of the perpendicular magnetic energy spec- 
trum determines the energy spectra of the other fields. 
This linearity assumption appears to be well satisfied 
in comparisons to nonlinear kinetic simulation results 8 . 
This approach enables us to fit three different curves by 
adjusting only the two Kolmogorov constants, C\ and 
C2, in the weakened cascade model, providing increased 
confidence in fits to numerical spectra. 
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FIG. 3. Energy spectra from the Pi — 1 dissipation range 
simulation^ resolving both the ion and electron Larmor ra- 
dius scales, depicted by vertical dotted lines. For a pi = 1 
and Ti/T c = 1 plasma, thick lines present numerical energy 
spectra for the perpendicular magnetic (black solid), electric 
(green dashed), and parallel magnetic (magenta dot-dashed) 
fields. Predicted energy spectra from the nonlocal model (thin 
solid) for Ci = 1.65 and the local model (thin dashed) for 
Ci = 1.8 are overplotted for comparison. 



1. Moderately Damped ft = 1 Case 

The first comparison of the local and nonlocal weak- 
ened cascade models tests their ability to fit the results of 
the j3i = 1 dissipation range simulation, the first kinetic 
simulation resolving both the ion and electron Larmor ra- 
dius scales in a single simulation, the full details of which 
are reported in a companion paper—. The plasma param- 
eters for this gyrokinetic simulation using AstroGK— are 
Pi = 1 and Ti/T e — 1, and all dissipation is provided 
by physically resolved collisionless damping via the Lan- 
dau resonance onto the ions and electrons. In Figure [3] 
are plotted the one-dimensional energy spectra (thick 
lines) for the perpendicular magnetic field fluctuations 
Eb ± (black solid) , the parallel magnetic field fluctuations 
Eb« (magenta dot-dashed), and the perpendicular elec- 
tric field fluctuations E E± (green dashed). 

For comparison, the predicted energy spectra from the 
the nonlocal weakened cascade model (thin solid) and 
the local cascade model (thin dashed) are overplotted on 
Figure [3] Using the linear gyrokinetic eigenfunctions for 
the Alfven mode enables the determination of the paral- 
lel magnetic and perpendicular electric field spectra from 
the perpendicular magnetic spectrum output by the cas- 
cade model. For the nonlocal model, the Kolmogorov 
constants required to yield a good fit to the numerical 
simulation results are C\ = 1.65 ± 0.20 and C2 = 1.0, 
while for the local model they are C\ = 1.8 ± 0.35, 
C2 = 1.0, and C3 — 2.25. Note that a higher value of C\ 
leads to stronger weighting of the linear damping relative 
to the nonlinear energy transfer. The Kolmogorov con- 
stant C\ is the primary adjustable parameter in the weak- 
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FIG. 4. From the local and nonlocal weakened cascade models 
for the ft = 1 and Ti/T e = 1 plasma depicted in Figure O 
(a) the energy cascade rate e/eo vs. k±pi for local (dotted) 
and nonlocal (dashed) models, (b) the fractional contribution 
to e due to the large-scale shearing motions e a /e (blue), the 
local-scale motions e;/e (green), and the small-scale diffusive 
motions e^A (red) with linear kinetic damping (solid) and 
with no damping (dotted), and (c) the nonlinear frequency 
ui n i from the local model (dotted) and the nonlocal model 
(dashed). 



ened cascade model, dominantly controlling the shape 
of the energy spectrum. The second Kolmogorov con- 
stant Ci fine tunes the condition of critical balance, and 
has not been adjusted; tests of the distribution of en- 
ergy in wavevector space are necessary to constrain this 
Kolmogorov constant and are beyond the scope of the 
present work. 

Comparison of the energy spectra indicates that both 
the nonlocal and local cascade models are able to repro- 
duce the (3i — 1 simulation spectra with similar values 
for C\. Differences between the models become clear as 
we look more closely at various contributions to the en- 
ergy cascade rate, as presented in Figure |H In panel (a), 
the energy cascade rate e is plotted vs. k±pi for local 
(dotted) and nonlocal (dashed) models. This compari- 
son shows little difference between models, so we must 
look more closely at the local and nonlocal contributions 
to the energy cascade rate. 

In panel (b) of Figure HI for the nonlocal model, the 
fractions of the energy cascade rate from the large-scale 
shearing motions e s /e (blue), the local-scale motions e;/e 
(green), and the small-scale diffusive motions e^/e (red) 
are plotted. To highlight the effects of kinetic dissipa- 
tion on these contributions to the energy cascade rate, 
dotted lines give the results when kinetic dissipation is 
artificially set to zero. From this plot, it is clear that, 
as the cascade proceeds to higher wavenumber, the dif- 
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fusive contribution (red) to the cascade diminishes first 
at k^pi ~ 1 due to dissipation (compared to the un- 
damped case given by the red dotted line), leading to a 
fractional increase in the local and shear contributions. 
Next, the fraction of e due to local motions (green) begins 
to diminish at around k±pi ~ 20, leading eventually to 
a dominance of the energy cascade rate by the shearing 
motions of the large scales (blue). It is this dominance 
of the energy cascade rate by the large scale motions as 
kinetic damping dissipates the turbulent motions that is 
the primary difference between the local and nonlocal 
models. Note that the unusual peak at the left for e;/e is 
due to the window that defines the local contributions 61 . 

The difference between the local and nonlocal mod- 
els can be seen at high wavenumbers in the nonlin- 
ear frequency U) n i, plotted in panel (c). For the local 
model (dotted), the nonlinear frequency peaks at around 
k±Pi ~ 50, and then begins to diminish due to the dis- 
sipation of the local motions responsible for the nonlin- 
ear energy transfer. The nonlocal model (dashed), on 
the other hand, merely flattens out, as large scale mo- 
tions continue to support the nonlinear energy transfer 
at smaller scales. 

In summary, both the local and nonlocal models yield 
similar results in modeling the turbulent energy spectra 
in the moderately damped /3j = 1 shown in Fig- 

ure [3) The differences become apparent only as the ki- 
netic dissipation becomes sufficiently strong to diminish 
the local contribution to the nonlinear energy transfer, 
enabling nonlocal, large-scale shearing motions to domi- 
nate the nonlinear frequency, as seen at the high k±pi end 
of panel (b) in Figure [4] It is this difference in the phys- 
ical mechanisms that will prove crucial in cases where 
the kinetic damping is stronger, requiring the additional 
physics of the effect of nonlocal motions on the energy 
transfer to model correctly the steady-state energy spec- 
tra. 



2. Strongly Damped = 0.01 Case 

In a low beta plasma, the kinetic damping of fluctua- 
tions in the KAW regime is substantially stronger. In this 
more strongly damped case, the difference between the 
local and nonlocal models is dramatic: the local model 
is simply unable to fit the shape of the spectrum. In 
this section we compare the spectra predicted by the lo- 
cal and nonlocal cascade models to the steady state of 
the /3i = 0.01 nonlinear AstroGK simulation. In Fig- 
ure 03 panel (a) shows a fit of the nonlocal model to the 
AstroGK simulation spectra (same legend as Figure [3j, 
using Ci = 2.85 (thin solid) and with C\ = 2.85 ± 0.15 
(thin dashed). In panel (b) is presented the best fit us- 
ing the local cascade model, for C\ = 2.6 (thin solid) 
and C\ = 2.6 ± 0.4 (thin dashed). All cascade models 
in this figure use Ci = 1 and, as usual, the local model 
employs C3 = 2.25 so that the Kolmogorov constants 
of both models are comparable. It is immediately ap- 
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FIG. 5. Energy spectra from the nonlinear AstroGK gyroki- 
netic simulation of turbulence in a plasma with /3i = 0.01 and 
Ti/T e — 1. Thick lines represent numerical energy spectra 
for the perpendicular magnetic (black solid), electric (green 
dashed), and parallel magnetic (magenta dot-dashed) fields, 
(a) The nonlocal model with Ci = 2.85T0.15, where thin solid 
lines are the spectra for the central value, and the dashed lines 
demonstrate the ±0.15 variation, (b) The local model with 
Ci = 2.6 ± 0.4. 



parent that the local model is incapable of fitting the 
correct shape of the spectrum, generally showing power 
law slopes that are too flat at the low wavenumbers and a 
cutoff that is too sharp at high wavenumbers; the neglect 
of the effect of nonlocal motions on the energy cascade 
rate at high wavenumbers, where the kinetic dissipation 
becomes strong, leads to this failure of the local model. 

An inspection of the nonlinear frequency u) n i for both 
models, plotted in Figure [BJ further illustrates this point. 
The nonlinear frequency for the local model (dotted) 
peaks at about k±pi ~ 6, and then drops off rapidly. 
This occurs because strong kinetic damping dissipates 
the turbulent fluctuations at the local scale, consequently 
slowing the nonlinear energy transfer due to those local 
motions and enabling the linear kinetic damping to dom- 
inate over energy transfer at that scale, resulting in a 
sharp cutoff of the turbulent energy spectra. The non- 
local model (dashed), on the other hand, shows that the 
nonlinear frequency flattens to a constant value at high 
wavenumbers but does not decrease. In this case, it is the 
nonlocal, large-scale shearing motions 6 -, that dominate 
the nonlinear energy transfer rate at high wavenumbers, 
leading to turbulent spectra that do not cut off sharply 
and are able to fit the nonlinear numerical results. This 
evidence suggests that it is the effect of nonlocal motions 
on the nonlinear energy transfer rate that is responsible 
for the more slowly dissipating, nearly power-law appear- 
ance of the turbulent spectra in both recent solar wind 
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FIG. 6. The nonlinear frequency u n i for the local (dotted) 
and nonlocal (dashed) cascade models for the Pi = 0.01 and 
Ti/Te — 1 turbulent plasma shown in Figure [5] 

observations^— and nonlinear numerical simulations 9 . 

In summary, although the local model can fit the data 
in certain moderately damped cases, the effect of nonlo- 
cal, large-scale shearing motions on the nonlinear energy 
transfer deep in the dissipation range is essential to avoid 
a sharp cutoff of the spectra and fit the nearly power-law 
behavior observed in both solar wind observations and 
nonlinear numerical simulations 



3. Weak Dissipating KAW Turbulence (WDKT) 

A graphical depiction of the contribution of nonlocal 
turbulent fluctuations to the energy cascade rate at a 
particular wavenumber k±_ illuminates the effect of non- 
locality. Here the weakened cascade model, defined in 
EIIII CI is used to solve for the steady state of a plasma 
with fii = 1 and Ti/T e — 16, parameters chosen so 
that there is a sufficient dynamic range of the kinetic 
Alfvcn wave regime to realize the asymptotic limit of the 
undamped KAW cascade before reaching k±p e ~ 1 at 
k±Pi ~ 170. The model covers a range of wavenum- 
bers k±pi G [10~ 3 , 10 3 ], employs constants C\ — 4.5 and 
C*2 = 1, and is driven in critical balance with x(^-Lo) = 1. 
In Figure the perpendicular magnetic energy spec- 
trum Eb±_ in panel (a) reveals a strong MHD turbulence 
regime over k±pi 6 [10~ 3 , 1] and a strong KAW turbu- 
lence regime over k±pi G [1, 10 2 ]. The kinetic dissipation 
begins to significantly alter the spectra for k±pi > 100, 
leading to a transition from strong to weak turbulence 
locally. But this situation is not a standard weak tur- 
bulence picture, such as those typically considered in 
incompressible MHD turbulence. Rather, the presence 
of large-scale shearing motions and the importance of 
further kinetic dissipation leads to a very different kind 
of weak KAW turbulence: we refer to this as the range 
of weak dissipating KAW turbulence (WDKT), seen over 
wavenumbers k±pi G [10 2 ,10 3 ] in panel (a) of Figure[7] 

Plotted in the lower three panels of Figure [7] is the 
function in the integrand of equation (|26|) for the total 
nonlinear frequency uj n i(k±) at (b) k±pi = 0.0316, (c) 
k±pi = 10, and (d) k±pi = 200, where the value of local 
wavenumber k± is indicated in each plot by the verti- 




FIG. 7. (a) The perpendicular magnetic energy spectrum for 
strong turbulence in a plasma with Pi = 1 and Ti/T e — 16 
predicted by the weakened cascade model. In the three 
lower panels is plotted the function in the integrand of equa- 
tion (|26[) for the total nonlinear frequency u) n i(k±) at (b) 
k ± pt = 0.0316, (c) k ± pi = 10, and (d) k ±Pi = 200. The 
value of the local wavenumber k± is indicated in each plot by 
the vertical solid line, and shading represents contributions 
due to motions on local scales from k±/2 to 2k±. 



cal solid line. The contributions to the integral yielding 
Wni(fcx) due ^0 local motions is given by the shaded re- 
gion below the curve (from k±/2 to 2k±); below the curve 
at lower wavenumbers is the contribution from large- 
scale shearing motions, and at higher wavenumbers, from 
small-scale diffusive motions. In panel (b) is the con- 
tribution to the nonlinear frequency at k±pi = 0.0316 
in the strong MHD turbulence regime. The nonlinear 
frequency in this regime is dominated by local motions, 
with large-scale shearing contributing slightly more than 
small-scale diffusive motions. In panel (c) is the contri- 
bution to the nonlinear frequency at k±pi — 10 in the 
strong KAW turbulence regime. Again, u> n i(k±) is domi- 
nated by local motions, and here small-scale diffusive mo- 
tions contribute slightly more than large-scale shearing. 
In panel (d) is highlighted one of the key physics points 
of the weakened cascade model: the nonlinear frequency 
at k±pi = 200 in the weak dissipating KAW turbulence 
regime is dominated by large-scale shearing motions and 
not by local scale motions. Additionally, there is no con- 
tribution from small-scale diffusion because the cascade 
is terminated and no motions exist at smaller scales. 



C. Complete Spectrum 

The power of the weakened cascade model is best sum- 
marized by a single final example that demonstrates all 



14 



of the physics incorporated, including the transition from 
weak to strong MHD turbulence and the complementary 
transition from strong KAW to weak dissipating KAW 
turbulence. This case models the entire turbulent cas- 
cade over k±pt E [1CT 5 , 300] for a ft = 1 and T l jT e = 9 
plasma with weak energy injection %(fcxo) = 0.1 at 
k±oPi — k\\QPi — 10~ 5 . The Kolmogorov constants are 
Ci = 1.4 and C 2 = 1.0. 

The steady-state perpendicular magnetic energy spec- 
trum Eb ± for this case is presented in panel (a) of Fig- 
ure [H Panel (b) contains the evolution of the paral- 
lel wavenumber k\\pi and panel (c) shows the evolution 
of the nonlinearity parameter \ over the entire cascade. 
The solution shows weak MHD turbulence over the range 
k±pi s [10 -5 , 4 x 10 -4 ], with a one-dimensional magnetic 
energy spectrum Eb ± oc kj 2 , no parallel cascade, and 
weak nonlinearity x < 0-5- Then comes a regime of tran- 
sition from weak to strong MHD turbulence in the range 
k±p t £ [4 x 10~ 4 ,4 x 10" 3 ]. Over k±p t G [4 x 10- 3 ,1] 

— 5/3 

is strong MHD turbulence with spectrum Eb ± oc k ± , 
parallel cascade scaling as fcii oc k 2 ^ 3 according to crit- 
ical balance, and a nonlinearity parameter x ~ 1- The 
transition from strong MHD turbulence to strong KAW 
turbulence occurs at k±pi ~ 1. The range of strong 
KAW turbulence over k±pi e [1,30] yields a spectrum 

Eb ± oc kj 7 ^ 3 , parallel cascade scaling as fcii oc k 1 ^ 3 ac- 
cording to critical balance, and a nonlinearity parameter 
X ~ 1. Finally as the kinetic dissipation begins to be- 
come significant, the cascade becomes weak dissipating 
KAW turbulence (WDKT) over k ±P% e [30,300], with a 
spectrum that drops off exponentially, inhibition of the 
parallel cascade, and a dropping nonlinearity parameter 
X< 1. 



VI. DISCUSSION 

The weakened cascade model, summarized in §111 C[ is 
a refinement of an earlier cascade model^ intended to ex- 
plain better the spectra observed in nonlinear numerical 
simulations^ and recent high time resolution observations 
of the dissipation range of turbulence in the near Earth 
solar wind^~— . The new physical effects incorporated 
into this model are (1) the transition between weak and 
strong turbulence and (2) the effect on the nonlinear tur- 
bulent energy transfer by nonlocal fluctuations. 

The model is constructed to have as few free parame- 
ters as possible, letting the linear physics of the kinetic 
plasma dictate the character of the turbulent fluctua- 
tions. The weakened cascade model has only two free 
parameters in the form of dimensionless, order-unity Kol- 
mogorov constants: C\, which adjusts the weighting of 
the nonlinear energy transfer to the linear kinetic damp- 
ing, and C*2, which fine tunes the condition of critical 
balance in the nonlinearity parameter x- To determine 
a value of C2 from numerical simulations or observations 
is beyond the scope of the present work, so in all cases 



we have set C% = 1, essentially leaving only a single 
adjustable parameter C\ to determine the shape of the 
steady-state energy spectra. As presented in §V] this sin- 
gle degree of freedom in the weakened cascade model is 
sufficient to fit closely the shape of the turbulent spec- 
tra from nonlinear gyrokinetic simulations for both mod- 
erately damped ft = 1 and strongly damped ft = 0.01 
cases, giving us confidence that the model contains the es- 
sential ingredients necessary to describe successfully the 
energetics of the turbulent cascade in a weakly collisional 
plasma. 

Two novel ingredients, inspired by turbulence phe- 
nomenology, are central to the weakened cascade model. 
The first is the prescription for the evolution of the paral- 
lel cascade depending on the strength of the nonlinearity 
parameter, given by equation (|3T))) . This equation is con- 
structed to reproduce the predictions of strong and weak 
turbulence theories in the appropriate limits. The sec- 
ond, and likely more important, ingredient is the nonlocal 
form of the nonlinear energy transfer frequency, given by 
equation (|2l)|) . The abandonment of the Kolmogorov hy- 
pothesis of locality^ is the advance necessary to explain 
the nearly power-law spectra seen in both the nonlinear 
simulation results^ and the solar wind dissipation range 
observations 1 ^— . When dissipation begins to weaken 
the local scale motions, the smaller scale motions are nec- 
essarily yet more strongly damped, so diffusion by those 
smaller scale motions is effectively negligible. Therefore, 
the key result of the weakened cascade model presented 
here is that the nonlocal effect of shearing by large scale 
motions explains the nearly power law appearance of nu- 
merical and observational spectra in the dissipation range 
(this effect is particularly critical in the low beta case pre- 
sented in WB2|) . 

The locality of MHD turbulence in the inertial range 
has gained a significant amount of attention in the lit- 
erature recently, with some studies finding evidence for 
nonlocalit y 48 i 63 ~ — , while other studies claim locality in 
the asymptotic limit of large Reynold's numbers 66 (sec 
Mininni 2010^ for a recent review of this topic). At 
question is whether the properties of the MHD inertial 
range become independent of the details of the mecha- 
nisms responsible for driving and dissipating the turbu- 
lence, enabling the determination of a self-similar solu- 
tion for a universal MHD turbulent spectrum. When the 
dissipation range of plasma turbulence is being consid- 
ered, however, it is obvious that no self-similar solution 
can be found, so the importance of a nonlocal contribu- 
tion to the energy cascade rate in the dissipation range, 
as proposed here, should come as no surprise. 

It is interesting to note that scale locality would be re- 
covered, and therefore a self-similar solution would arise, 
if the kinetic Alfven wave range of scales, pi <C 1/fc <C p e , 
becomes asymptotically large, as occurs in the limit of 
vanishing electron mass, m e — > 0. In this very small 
electron mass limit m e —¥ 0, the kinetic Alfven wave par- 
allel phase velocity is always much less than the electron 
thermal velocity, ui/ku <C v te , so collisionless damping via 
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FIG. 8. Numerical solution of the weakened cascade model for the steady state spectrum in a /3; = 1 and Ti/T e = 9 plasma. 
The panels plotted are the same as Figure [2] Regimes indicated in panel (a) are weak MHD turbulence, a transition regime 
from weak to strong MHD turbulence, strong MHD turbulence, strong KAW turbulence, and finally weak dissipating KAW 
turbulence (WDKT). 



resonant wave-particle interactions with the electrons be- 
comes negligible. In the absence of damping, a second in- 
ertial range is recovered for the dispersive kinetic Alfven 
waves over the range of scales pi <C 1/k <C p e - The weak- 
ened cascade model describes this limit correctly, because 
if m e — > 0, then the collisionless damping rate 7 — > in 
equation (fl~3|) . In this case, the weakened cascade model 
recovers a second inertial range of kinetic Alfven waves 



with a spectral index of —7/3, as shown by the dashed 
line denoted the "undamped model" in Figure 3 of the 
original cascade model paper by Howes et al. (2008) 5 . 
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A. Comparison to Previous Cascade Models 

It is instructive to compare the predictions of the weak- 
ened cascade model, with its inclusion of nonlocal effects 
on the energy cascade rate, to those of the original cas- 
cade model^ and of a similar model by Podesta et al. 
2010^, both strictly local models. Two of the general 
conclusions of the original cascade model were that vari- 
ations in collisionless damping with plasma parameters 
naturally explained observed variations in solar wind tur- 
bulent spectra and that the dissipation range spectrum 
should be an exponential fall ofl&. Podesta et al. make the 
stronger claim that "an energy cascade consisting solely 
of KAWs cannot reach scales of order the electron gyrora- 
dius, k±p e ~ 1. It implies that the power-law spectrum 
in the regime of electron scales must be supported by 
wave modes other than the KAW" 68 . 

Since the only wave mode supported in the dissipation 
range of a fii = 1 plasma in gyrokinetic theory is the 
kinetic Alfven wave, the = 1 gyrokinetic simulation 9 
demonstrates unequivocally by counterexample that the 
claim by Podesta et al^ — that the kinetic Alfven wave 
cascade cannot reach scales k±p e ~ 1 — is incorrect. We 
explain below that the primary reason for this failure lies 
in the determination of the Kolmogorov constant in their 
cascade model. 

Each of these cascade models incorporates an order- 
unity dimensionless constant that essentially adjusts the 
weighting of nonlinear energy transfer to that of the lin- 
ear kinetic damping. Podesta et al. estimate k±pi ~ 28 as 
the wavenumber at which the KAW cascade terminates, 
given by their equation (23), which is linearly dependent 
on their constant A. Without any guidance to determine 
this constant, Podesta et al. base their conclusions on the 
theory with A = 1. We can usefully compare their model 
with the local version of the weakened cascade model 
(using equation [3T] for the nonlinear frequency) to de- 
termine what value of A would fit the /3j = 1 dissipation 
range simulation spectra in Figure |3l Using b\ — kE(k) 
in the limit k^ 3> fe|| to connect to their model, the en- 
ergy cascade rate e in our local model is given by 

e = C; 3/2 C 3 J^ C) kE(k). (32) 

The energy cascade rate e in Podesta et al., based on 
their equation (6)^, is given by 

A 

e=—ujkE(k). (33) 

The assumption of critical balance is taken by Podesta 
et al. to be w = W £° C \ so the two models should give 
similar results if A = 2ttC\ C3. Substituting in the 
best fit values C\ — 1.8 and C3 = 2.25 from Figure El we 
find A = 5.85. This comparison suggests that Podesta 
et al. significantly underestimated the weight of the en- 
ergy cascade rate relative to the linear kinetic damping, 
leading to a conclusion inconsistent with direct numerical 



simulations of kinetic turbulence. Nonlinear kinetic sim- 
ulations of turbulence therefore play an invaluable role 
in the effort to understand turbulence in the solar wind 
and other weakly collisional astrophysical plasmas. 



B. Limitations of the Weakened Cascade Model 

The weakened cascade model has been developed to 
predict the turbulent energy spectra occurring in weakly 
collisional astrophysical plasmas. We discuss here a num- 
ber of assumptions that have been made in the construc- 
tion of the model. 

First, the model has been constructed to reproduce the 
scaling of the turbulent energy spectra given by the Gol- 
dreich and Sridhar theories for weak and strong MHD 
turbulenc o 24 i 25 i 31 , and their extension to turbulence in 
the kinetic Alfven wave regime. An alternative theory 
for the scaling of strong MHD turbulence has been pro- 
posed by Boldyrev^, and substantial numerical support 
has accumulated in its favo r 32 i 70 - ~— . Modifications of 
the weakened cascade model to reproduce instead the 
Bolydrev scalings will be discussed in a subsequent pa- 
per. 

Another potential limitation of the weakened cascade 
model is the conjecture that the parallel cascade of energy 
is inhibited for weak turbulence in both the MHD and 
KAW regimes, as discussed in §111 Al 

The one-dimensional nature of the weakened cascade 
model restricts its direct applicability to plasmas in which 
there is a single scale of energy injection at the outer 
scale (fc_Lo,fc||o)- As discussed in Howes et al. 2008^, the 
lack of structure at small scales in the solar wind energy 
spectrum is evidence against significant energy injection 
as scales smaller than the outer scale (fcj_Oj^||o)j so the- 
model appears to be broadly applicable to the solar wind. 

A number of other factors that may significantly af- 
fect the turbulent fluctuations in the solar wind are not 
incorporated into the weakened cascade model, includ- 
ing the radial expansion of the solar wind, kinetic tem- 
perature anisotropy instabilities, and the imbalance of 
sunward and anti-sunward Alfven wave energy fluxes. A 
detailed discussion of these effects is presented in the pa- 
per describing the original cascade model^ and will not 
be repeated here. It suffices to say that the weakened 
cascade model is an attempt to understand the quantita- 
tive details of the energy transport in balanced, Alfvenic 
turbulence in a weakly collisional plasma; the relation 
of these additional effects to this fundamental turbulent 
evolution merits further investigation. 

The weakened cascade model does not account for one 
recently discovered physical mechanism that may play 
an important role in the energy transport in weakly col- 
lisional Alfvenic turbulence: the entropy cascade£i22r— . 
As described in Schekochihin et al. 2009^, the ion en- 
tropy cascade is a dual cascade to small scales in both 
physical space and velocity space of the ion distribution 
function. Operating at scales below the ion Larmor ra- 
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dius k±pi > 1, in the regime of kinetic Alfven wave tur- 
bulence, this process is driven by nonlinear phase mixing 
and represents an alternative channel of energy trans- 
port that is not included in the weakened cascade model. 
Comparisons of the predictions of the weakened cascade 
model with the results of a suite of nonlinear gyrokinetic 
simulations will enable an evaluation of the importance 
of the ion entropy cascade in the turbulent energy trans- 
port, an important line of future research. 



VII. CONCLUSIONS 

Early cascade models for kinetic turbulence in weakly 
collisional astrophysical plasmas, such as the solar wind, 
suggested that the energy spectra in the dissipation range 
should exponentially fall ofS, and that due to collisionless 
damping, kinetic Alfven waves could not be responsible 
for the cascade of energy to electron scales^. However, 
nonlinear kinetic simulations of turbulence over the en- 
tire dissipation ranged and high time resolution obser- 
vations of the dissipation range fluctuations in the solar 
wind 10_ ~ yield nearly power-law energy spectra rather 
than an exponential decay. The failure of these early 
cascade models motivated refinements to explain the nu- 
merical and observational results, yielding the weakened 
cascade model presented here. 

The original cascade model by Howes et al. 2008 5 was 
based on three assumptions: (1) the Kolmogorov hypoth- 
esis of the locality of the nonlinear energy transfer in 
wavenumber space, (2) the conjecture of critical balance 
at all scales, and (3) the applicability of linear kinetic 
damping rates. The weakened cascade model eliminates 
the first two assumptions of the original model, result- 
ing in a more broadly applicable and physically realistic 
model. 

The assumption of critical balance at all scales is 
dropped; instead, the model handles explicitly the transi- 
tion between weak and strong turbulence. A key point in 
the treatment of weak turbulence is our conjecture that 
the parallel cascade is inhibited in both the MHD and 
kinetic Alfven wave regimes. This physics is contained in 
the model for the evolution of the parallel wavenumber 
given by equation (|30[) . The more significant advance 
of the weakened cascade model is the abandonment of 
the locality hypothesis of Kolmogorov. Although energy 
is still transferred locally in wavenumber space — for ex- 
ample, from k± to 2k± — the turbulent fluctuations re- 
sponsible for that energy transfer may be nonlocal. Both 
shearing by motions at larger scales and diffusion by mo- 
tions at smaller scales contribute to the nonlinear energy 
cascade rate according to equation (|26[) . 

In fjvl we have demonstrated that the weakened cas- 
cade model reproduces the transition from weak to strong 
MHD turbulence as predicted by theory. As the colli- 
sionless dissipation in the kinetic Alfven wave regime be- 
comes significant, the model also shows a complementary 
transition from strong KAW turbulence to weak dissipat- 



ing KA W turbulence, a new regime of weak turbulence in 
which the effect of shearing by large scale motions and 
continued kinetic dissipation play an important role. 

The key result of this paper is that the nearly power- 
law energy spectra observed in the dissipation range of 
both numerical simulations and solar wind observations 
are explained by the inclusion of the effect of nonlocal 
motions on the nonlinear energy cascade rate, specifi- 
cally the shearing by large-scale motions. Although nu- 
merical spectra for a moderately damped fii = 1 plasma 
may be equally well explained by either local or non- 
local models (Figure [3]) , for the more strongly damped 
Pi = 0.01 plasma, the inclusion of nonlocal effects is crit- 
ical for the model to fit the numerical energy spectra (Fig- 
ure [5]) . The importance of the nonlocal shearing motions 
to the energy cascade rate is demonstrated in panel (d) of 
Figure where it is clear that the large-scale contribu- 
tion dominates over the local contribution. The effect at 
these strongly dissipative scales is that the nonlinear fre- 
quency does not decrease with increasing perpendicular 
wavenumber, as a local model would suggest, but that it 
remains constant due to the large-scale contribution, as 
shown in Figure[6l Thus, by abandoning the Kolmogorov 
hypothesis of locality, the weakened cascade model ex- 
plains the nearly power-law spectra found in numerical 
and observational studies of the dissipation range by in- 
cluding the nonlocal effect of large-scale shearing motions 
on the energy transfer rate. 

The ultimate aim of the weakened cascade model is 
not to fit numerical and observational turbulent spectra, 
but to predict them. Nonlinear kinetic simulations of 
dissipation range turbulence have already played an im- 
portant role in evaluating the weakened cascade model. 
Since the Kolmogorov constants C± and C2 may, in prin- 
ciple, depend on the plasma parameters fti and Ti/T e , 
numerical studies will continue to play a critical role as 
we perform a suite of kinetic turbulence simulations over 
a range of plasma parameters to determine this depen- 
dence, ultimately striving for predictive capability. In 
addition, these numerical studies will also be crucial in 
judging the importance of the ion entropy cascade to the 
turbulent energy transfer in the dissipation range. 
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Appendix A: The Effect of Nonlocal Interactions on 
Turbulence Scaling 

Here we consider the effect of the nonlocal contribution 
to the nonlinear frequency, given by equation (f26|) . on 
the scaling of the spectra in the strong MHD and strong 
KAW inertial ranges. 

In the strong MHD inertial range, the contribution to 
the nonlinear frequency due to motions at each wavenum- 

ber is w^ ; 0C ' ) = k±bk = Cj^e 1 / 3 /^ 3 . Performing the 
integral in equation (|26p to find u> n i, we obtain 



-Lmax } ' 



(Al) 



where the first term is due to the large-scale shearing, and 
the second is due to the small-scale diffusion. The inertial 
range is defined as the range of scales unaffected by large- 
scale driving or small-scale dissipation, corresponding to 
the limit k±Q <C fcj_ <C k± max . In this limit, the nonlinear 
frequency simplifies to 



Unl(kl 



3 (loc) 3 (be) _ 9 (loc) 



-a 



(A:j 



(A2) 



Therefore, within the strong MHD inertial range, the 
nonlinear model yields a nonlinear frequency of the form 

w n i(k±) — CsUJni >C \k±), where the constant is C3 = 
2.25. This order-unity constant factor is the only differ- 
ence between models, so the scaling of the nonlocal model 
within the strong MHD inertial range will therefore be 
the same as a local model. Note that the contribution 
of the large-scale shearing motions is twice that of the 
small-scale diffusive motions. 

In the strong KAW inertial range, we have w^ 02 ' = 

kj_bkk±pi = 'e 1 / 3 p^ 3 kV 3 . Performing the integral 
in equation (1261) . we find 

^)=\ci^' 3 P r{kT-k^ 
+ \c\'\^ P rki{k^ 3 -kii{i)^) 

where the first term is the large-scale contribution and 
the second is the small-scale contribution. Within the 
KAW inertial range, we apply the limit k± <C k± <C 
k±max to simplify the nonlinear frequency to 



3 (loc) 3 (loc) _ 9 (loc) 



OJnl{k±) = ~UJ nl +2 W rU 



(A4) 



Again, we find that the scaling of the nonlocal model in 
the strong KAW inertial range will be that same as that 
of a local models, the only difference being the same con- 
stant factor C3 = 2.25. In this case, however, the small- 
scale diffusive motions contribute twice as much to the 
nonlinear frequency as the large-scale shearing motions. 
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